Suppression of Giant Magnetoresistance by a superconducting contact 
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We predict that current perpendicular to the plane (CPP) giant magnetoresistance (GMR) in 
a phase-coherent magnetic multilayer is suppressed when one of the contacts is superconducting. 
This is a consequence of a superconductivity-induced magneto-resistive (SMR) effect, whereby the 
conductance of the ferromagnetically aligned state is drastically reduced by superconductivity. To 
demonstrate this effect, we compute the GMR ratio of clean (Cu/Co)„Cu and (Cu/Co) n Pb multi- 
layers, described by an ab-initio spd tight binding Hamiltonian. By analyzing a simpler model with 
two orbitals per site, we also show that the suppression survives in the presence of elastic scattering 
by impurities. 



PACS: 75.50Pa, 74.80Dm, 72.00 

During the past decade, electronic properties of hy- 
brid nanostructures have been widely studied, both from 
a fundamental point of view and for their potential ap- 
plications. At a fundamental level, new physics associ- 
ated with such structures arises from the proximity of two 
electronic ground states with different correlations, which 
can reveal novel scattering processes not apparent in the 
separate materials. For example normal-superconducting 
hybrids, formed when a normal (N) metal or semiconduc- 
tor is placed in contact with a superconductor (S), have 
been shown to exhibit a range of unique transport phe- 
nomena 0J^] associated with Andreev scattering at the 
N-S boundary. Similarly ferromagnetic (F)-normal mul- 
tilayers and spin-valves exhibit magnetoresistance prop- 
erties HJiJ] associated with spin-filtering by the F-layers. 
Perhaps the most spectacular of these is giant magne- 
toresistance (GMR) in magnetic multilayers, whose re- 
sistance, with the current perpendicular to the planes of 
the layers (CPP), can change by more than 100% un- 
der the application of a magnetic field. Until recently 
N-S and F-S nanostructures have been studied in isola- 
tion, but during recent months, a number of experiments 
have demonstrated that N-F-S hybrids exhibit a range 
of novel features |^-^|, including a long-range supercon- 
ducting proximity effect in the F- material, extending over 



length scales for an excess of the magnetic length y jP~ , 

where D is the diffusion coefficient and E ex the exchange 
splitting. 

The aim of this Letter is to address a new phenomenon, 
namely the effect of superconductivity on CPP-GMR in 
phase-coherent multilayers of the type (N/F)„S, where S 
is either a superconductor or a normal metal. We shall 
demonstrate that as superconductivity is induced in the 
S-contact (eg by lowering the temperature) CPP-GMR 
is almost completely suppressed. This result is remark- 
able, since for example CPP-GMR experiments by the 
Michigan State University group already employ 

superconducting contacts. We shall argue that GMR in 
such experiments is a consequence of spin-flip scattering, 
which if eliminated, would cause a dramatic suppression 
of GMR. To demonstrate the superconductivity-induced 
suppression of GMR, we use the method outlined in JO] 
to compute the zero-bias, zero-temperature conductance 
of the (Cu/Co)„Pb multilayer sketched in figure 0, de- 
scribed by an spd tight-binding Hamiltonian, with tight- 
binding parameters fitted to accurate ab-initio density 
functional calculations jl^ ]. Since we are interested in a 
phase-coherent nanostructure in which the magnetic mo- 
ments of successive Co layers are either parallel or anti- 
parallel, the conductance in the normal state is given by 
the Landauer formula 
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where T a — Tr t a t a \ with t the multi-channel trans- 
mission matrix for the structure. In the superconducting 
state, equation ( jl) is replaced by current-voltage rela- 
tions derived in and re-derived in fll4] , [i"5| |, which in 
the absence of quasi-particle transmission through the 
superconductor yields 



Gns 



4e 2 



-Ra 



(2) 



where R a = r°r°* is the Andreev reflection coeffi- 
cient, which for a spin-singlet superconductor is inde- 
pendent of the spin a of the incident quasi-particle. In 
what follows, when the magnetic moments of adjacent 
Co layers are aligned (anti-aligned) we denote the con- 

and therefore the 



ductances by G£j N , G£j S 



(G NN , G NS ) 



GMR ratios are given by Mnn 
Mns = (Gns - Gns)/ g ns- 
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Consider hrst the case of quasi-ballistic transport in 
which there is no disorder within the layers, nor at the 
interface, but the widths of the Co layers are allowed 
to fluctuate randomly by 1 atomic layer. Such a struc- 
ture is translationaly invariant in the direction par- 
allel to the layers, and the Hamiltonian is diagonal in 
a Bloch basis (fcii). Therefore the trace over scattering 
channels in equation (|f|) and (|J) can be evaluated by 
computing the scattering matrix at separate fcii points. 
Figure 2a shows results for the GMR ratio in the normal 
and superconducting states, obtained by summing over 
5 • 10 3 fen points, and clearly demonstrates a dramatic 
superconductivity-induced suppression of GMR. Figure 
2b and 2c show results for the individual conductances 
(note the difference in scales) and demonstrate that the 
GMR ratio Mns is suppressed because G^ s is drastically 
reduced compared with G^ N . 

To understand this effect, consider the simplest possi- 
ble model of spin-dependent boundary scattering shown 
in figure [|, which in the limit of delta-function F lay- 
ers, reduces to the model used to describe the N-F- 
S experiment of ||. Fig 3a (3b) shows a cartoon of 
a majority (minority) spin, scattering from a series of 
potential barriers in successive aligned F layers. Since 
the minority spins see the higher barrier, one expects 
Tf < TF . Figures 3c and 3d show the scattering poten- 
tials for anti-ferromagnetically aligned layers, for which 
T A = T A < Tf. For such an ideal structure, GMR arises 
from the fact that TF ^> Tf and T^. In the presence of a 
single superconducting contact this picture is drastically 
changed. For ferromagnetically aligned layers, figure 3e 
shows an incident majority electron scattering from a se- 
ries of low barriers, which Andreev reflects as a minority 
hole and then scatters from a series of high barriers (fig- 
ure 3f). The reverse process occurs for an incident minor- 
ity electron, illustrating the rigorous result that the An- 



dreev reflection coefficient is spin-independent. Figures 
3g and 3h illustrate Andreev reflection in the anti-aligned 
state. The crucial point illustrated by these sketches is 
that in presence of a S contact for both the aligned (fig- 
ures 3e and 3f) and anti-aligned (figures 3g and 3h) states 
the quasi-particle scatters from N (=4 in the figures) high 
barriers and N (=4) low barriers and therefore at the level 
of a classical resistor model, one expects G^ s ss G^ s . 

Of course the rigorous results of figure ||, obtained us- 
ing an spd Hamiltonian with 36 orbitals per atomic site 
(spdx2 for spin x2 for particle- hole degrees of freedom) 
go far beyond this heuristic argument. In the case of 
aligned or anti-aligned F-layers the problem involves two 
independent spin fluids and therefore the Hamiltonian is 
block-diagonal with 18 orbitals per site. The Hamilto- 
nian used to obtain these results is of the form 
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where (Hr) describes a semi-infinite left-hand 
(right-hand) crystalline lead, #lm (#mr) is the coupling 
matrix between surface orbitals on the left (right) lead 
and the left (right) surface of the magnetic multilayer 
and Hm is the tight-binding Hamiltonian describing the 
multilayer. Consider first the retarded Green's function 
g = (E — Hl — Hr + iQ + )~ 1 of the two decoupled semi- 
infinite leads. If the surfaces of the leads each contain 
M atoms, then TJlm and Hmr are 36M x 36M matri- 
ces and the portion of g involving only matrix elements 
between orbitals on the left and right lead surfaces is a 
(2 x 36M) x (2 x 36M ) block diagonal matrix g s whose 
matrix element gfj vanish for i,j belonging to different 
leads. Using a semi-analytic form for g s derived in [p"i"|| , 
the surface Green function G for the leads plus multi- 
layer can be computed by first recursively decimating the 
Hamiltonian H LM + H M + H M r to yield a (72M x 72M) 
matrix of couplings Hm between surface orbitals of the 
leads, and then computing the inverse 



G b = [(g*)- 1 - H M ] 
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Once the full surface Green's function G s is known, 
the scattering matrix elements between open scattering 
channels are obtained using a generalized Fisher-Lee re- 
lations Pq] . For the calculation of figure |^, involving fee 
crystalline leads aligned along the (110) direction the spd 
Hamiltonians for the bulk materials are known ]l2[ ] , but 
the hopping parameters at the interfaces between differ- 
ent materials are not currently available. As a simplest 
approximation, these surface couplings were chosen to be 
the geometric mean of their bulk values. 

Despite our use of a highly efficient recursive Green's 
function technique to exactly evaluate the scattering ma- 
trix of a multilayer, currently available computing re- 
sources restrict such a calculation to systems with trans- 
lational invariance parallel to the planes. To demonstrate 
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that the suppression of CPP-GMR is a generic feature of 
N-F-S hybrids and to study the effect of elastic impurity 
scattering, we now examine a reduced two band (s-d) 
model with a Hamiltonian matrix 
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In this model, h 



a/3 



■ -i-JijSapSad with hi the exchange 
splitting on site i for the d orbital, which vanishes if i be- 
longs to a N or S layer and is of magnitude h if i belongs 
to a F layer; A"- = Ai5ij5 a /3 where the superconduct- 
ing order parameter Aj vanishes if i belongs to a N or F 
layer and equals A if i belongs to the S region (s-wave 
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otherwise. Note 



superconductivity); (H ) 

i,j nearest neighbors and (H a 
that this model is the minimal model including the pos- 
sibility of s-d interband scattering, that has been shown 
to play a crucial role in describing the scattering prop- 
erties of a transition metal multilayer [[nj. ef is chosen 
to be a random number, uniformly distributed between 
e a — w/2 and e Q + w/2. We choose the parameters of the 
model to fit the conductance and the GMR obtained from 



the spd model for Cu/Co 11 1, namely (all quantities are 



expressed in eV) e 



7& = 



0.85, 7& d u 
2-7, 7 C o 



Cu 



= -7.8 
1.1, e 



c Cu 



-4-0, 7& s u 



Co 



-4.6, e 



Co 



0.85, 7 g, = 0.9, h = 1.6, A 



= 2.7, 
-2.0, 
10" 3 , 



7co 
w = 0.6. 

Figure 4a shows results for the GMR ratios Mnn and 
A/ns and demonstrates that the suppression of CPP- 
GMR by superconductivity survives in the presence of 
disorder. We have investigated a range of higher disor- 
ders and system sizes and find superconductivity-induced 
GMR suppression in all cases. The disorder used in fig- 
ure H has been chosen to illustrate an additional novel 
feature, not so-far discussed in the literature, namely 
that ballistic majority spins can co-exist with diffusive 
minority spins. For a strictly ballistic structure, the con- 
ductance G is almost independent of length L and the 
product G ■ L varies linearly with L. This behavior oc- 
curs for the majority spin in the N-results of figure 4b. 
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In contrast, for minority spins the product G 
hibits a plateau for L > 500 AP, which is characteristic 
of diffusive behavior. The same plateaus are also present 
in the product G^ N • L and in the curves of G^s • L and 
G^ s • L shown in figure 4c. 



In summary we predict that the presence of a sin- 
gle superconducting contact destroys the sub-gap CPP- 
GMR of phase-coherent magnetic multilayer. This arises 
because superconductivity suppresses transport in the 
majority sub-band in the ferromagnetic alignment, but 
causes little effect in the antiferromagnetically aligned 
state. This drastic reduction in G^g compared with G^j N 
is itself a remarkable superconductivity-induced magneto 
resistance effect. This suppression will be lifted at high 
biases and finite temperatures, where transport occurs 
via both Andreev reflection and quasi-particle transmis- 
sion. Furthermore the presence of spin-flip scattering at 
the superconducting interface will destroy this effect, be- 
cause if the spin of an Andreev reflected hole is flipped 
by such process, before it traverses the multilayer, only 
the contribution to GMR from layers within a spin-flip 
scattering length of the N/S interface is suppressed. 
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FIG. 1. Scheme of the system considered: N-F multi- 
layer sandwiched between a normal (N) contact and a super- 
conducting (S) contact. In the calculation presented N=Cu, 
F=Co and S=Pb. 
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FIG. 2. GMR (a) and spin conductances of Cu/Co multi- 
layers as a function of the Cu thickness, in the case of normal 
contacts (b) and with one N contact replaced by a supercon- 
ducting (c) contact. The system is disorder-free in the layers, 
but the Co thickness is allowed to randomly fluctuate by 1 
atomic plane. Every point corresponds to the mean of 10 
different configurations. The error bars denote the standard 
deviation of the mean. 
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FIG. 3. Cartoon of the different scattering processes. Fig- 
ures 3a, 3b, 3c and 3d describes the transmission of spin elec- 
trons e^' in a NN system. Figures 3e, 3f, 3g and 3h describe 
the NS case. Note that in the F case a majority (minority) 
spin electron e 1 (er) is Andreev reflected as a minority (ma- 
jority) hole h. (h In the anti-aligned (A) case the path of 
the incoming electrons and outcoming holes is identical for 
both spins. The total number of large barriers is the same in 
the A and F case, and this produces GMR suppression. 
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FIG. 4. GMR (a) and spin conductances of Cu/Co mul- 
tilayers calculated with the s-d two bands model, in the case 
of normal contacts (b) and with one N contact replaced by 
a superconducting (c) contact. The Co and Cu thicknesses 
are fixed and are respectively 15 and 8 atomic planes. Every 
point on the graph corresponds to an additional double bi- 
layer Co/Cu/Co/Cu. The on-site energy fluctuates randomly 
according with a normal distribution of width w = 0.6, and 
the error bars are the standard deviation of the mean over 10 
random configurations. The unit cell is a square with 5x5 
atomic sites, and we consider 25 k— points in the Brillouin 
zone. The horizontal line denotes GMR=0. 
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